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We study mesoscopic transport in the quasi-one-dimensional wires and rings made of a two-dimensional 
conductor of width W and length L ^ W . Our aim is to compare an impurity-free conductor with grain 
boundaries with a grain-free conductor with impurity disorder. A single grain boundary is modeled as a set of 
the two-dimensional-(5-function-like barriers positioned equidistantly on a straight line and disorder is emulated 
by a large number of such straight lines, intersecting the conductor with random orientation in random posi- 
tions. The impurity disorder is modeled by the two-dimensional (5-barriers with the randomly chosen positions 
and signs. The electron transmission through the wires is calculated by the scattering-matrix method, and the 
Landauer conductance is obtained. Moreover, we calculate the persistent current in the rings threaded by mag- 
netic flux: We incorporate into the scattering-matrix method the flux-dependent cyclic boundary conditions and 
we introduce a trick allowing to study the persistent currents in rings of almost realistic size. We mainly focus 
on the numerical results for L much larger than the electron mean-free path, when the transport is diffusive. If 
the grain boundaries are weakly reflecting, the systems with grain boundaries show the same (mean) conduc- 
tance and the same (typical) persistent current as the systems with impurities, and the results also agree with the 
single-particle theories treating disorder as a white-noise-like potential. If the grain boundaries are strongly re- 
flecting, the rings with the grain boundaries show the typical persistent currents which can be about three times 
larger than the results of the white-noise-based theory, thus resembling the experimental results of Jariwala et 
al., Phys. Rev. Lett. 86, 1594 (2001). Finally, we extend our study to the three-dimensional conductors with 
columnar grains. We find that the persistent current exceeds the white-noise-based result by another one order 
of magnitude, similarly as in the experiment of Chandrasekhar et al., Phys. Rev. Lett. 67, 3578 (1991). 

PACS numbers: 73.23.-b, 73.23.Ra 



I. INTRODUCTION 

Magnetic flux $ piercing the opening of a mesoscopic con- 
ducting ring gives rise to the equilibrium electron current cir- 
culating along the ring. This current is known as persistent 
current |1]. A single electron at the energy level En carries 
the current /„ = - dEn{^)/d^ 121,131]. At zero temperature 
the persistent current in the ring is given as / = ^ /„, where 
one sums over all occupied states below the Fermi level iQ]. 
In a single-channel ballistic ring, the amplitude of the persis- 
tent cuiTent is /o = evp/L, where e is the electron charge, vp 
is the electron Fermi velocity and L is the ring circumference. 
The current changes its sign when a single electron is added 
into the ring. In the ballistic ring with Nc conducting channels 
the amplitude of the persistent current scales as ^/T^Iq due to 
the random sign of the current in each channel |5]. 

If disorder is present, the persistent current does not vanish 
and the amplitude and sign depend on the specific configura- 
tion of disorder. Therefore, it is customary to study the typ- 
ical persistent cuiTent Ityp — (J^)^/^, where (. . . ) is the en- 
semble average. The authors of the works ^ analyzed the 
persistent cuiTent in a disordered normal-metal multi-channel 
ring by considering the non-interacting electron gas. They 
assumed that the electrons are scattered by disorder with the 
scattering potential ^(r) obeying the white-noise condition 
{V{r)V{r')) (X 6{r — r'). Using the Green's functions, they 
found for the typical current at $ = ± 0.25/i/e the formula 
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where the factor of 2 is due to the spin, / is the electron mean 
free path, ^ is the localization length, d is the dimensionality 
of the ring, and the condition Z <C i <C ^ means the diffusive 
regime. (The factor 1.6/d is derived in the appendix A from a 
more general formula from the literature.) 

The persistent current in a single isolated ring was for the 
first time measured by Chandrasekhar ef fl/. fs']. In that exper- 
iment, three different Au rings of size L ^ 100/ showed the 
persistent currents ranging from ^ 0.2ewi?/L to ^ 2evp/L, 
which is one-to-two orders more than predicts the formula 
Ifyp"^ ~ [evp / L){1/ L). This disagreement has sofar not 
been explained, specifically, the effort to explain it by con- 
sidering the electron-electron interaction (reviewed e.g. in 
|9]) was not successful. Ten years after the work |8], the 
same laboratory [10] prepared a new Au samples and ob- 
served the typical persistent currents much closer to the for- 
mula ~ {evF/L){l/L), but stiU two-to-three times 
larger. 

On the other hand, recent measurements of the persistent 
current in a single ring fTT], performed for thirty Au rings, 
have shown a good agreement with the formula ([T]). In ad- 
dition, the persistent current in a single Al ring has quite re- 
cently been measured by a new highly-sensitive method lll2ll . 
This work definitely demonstrates agreement of the experi- 
mental data with the formula ([T]l, modified by a temperature- 
dependent factor The experiments fil llll2ll thus strongly sug- 
gest, that the typical persistent current in a single disordered 
normal-metal ring is not affected by the electron-electron in- 
teraction at least for such metals like Au and Al. If this is the 
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case, then the disagreement between the formula ([T]l and pre- 
vious measurements of the Au rings is not due to the 
electron-electron interaction and the explanation, if any, may 
be hidden in the single-particle interaction with disorder 

The formula ([U holds for disorder modeled by the white- 
noise-like potential with spatially homogenous randomness. 
In reality, fabrication of the metallic wires/rings from such 
metals like Au, Ag, Cu, etc., involves techniques like the elec- 
tron beam lithography, lift-off, and metal evaporation, which 
provide wires/rings with disorder due the grain boundaries, 
impurity atoms and rough edges [iSJ. Thus it seems reason- 
able to study realistic disorder and to compare the results with 
the white-noise-based theory [6, 7]. 

In this work, electron transport in the mesoscopic wires and 
rings is studied with the aim to compare an impurity-free sys- 
tem with grain boundaries with a grain-free system containing 
the impurity disorder. (We ignore the edge roughness which 
is studied elsewhere ifo. [Til] .1 The mesoscopic wire is called 
quasi-one-dimensional (QID) if its length L is much larger 
than the width (W) and thickness (H) fil5il . We mainly study 
the QID wires and rings made of a two-dimensional conduc- 
tor (i/ 0) of width W and length L W, when the di- 
mensionality entering the formula ([T]l is d = 2 [7]. At the 
end we extend our study to the case d = 3, i.e., to the three- 
dimensional (3D) conductor with H ^ W. 

In our d = 2 study, a single grain boundary is modeled as a 
set of the two-dimensional-(5-function-like barriers positioned 
equidistantly on a straight line and disorder is emulated by a 
large number of such straight lines, intersecting the conductor 
with random orientation in random positions (figure [U. The 
impurity disorder is represented by many two-dimensional 
(5-barriers with randomly chosen positions and signs. The 
electron transmission through the wires is calculated by the 
scattering-matrix method lll6lll7ll . and the Landauer conduc- 
tance is obtained. To calculate the persistent current in the 
rings with magnetic flux, we include into the scattering-matrix 
method the flux-dependent cyclic boundary conditions and we 
introduce a trick allowing to study the typical persistent cur- 
rent in rings of almost realistic size. We mainly focus on the 
systems with L ^ I, when the transport is diffusive. 

If the grain boundaries are weakly reflecting, the systems 
with grain boundaries show for large enough L the same 
(mean) conductance and the same (typical) persistent current 
as the systems with impurities. The obtained results also agree 
with the single-particle theories 0, [Hi treating disorder as 
a white noise. 

If the grain boundaries are strongly reflecting, the rings with 
the grain boundaries are found to exhibit the typical persis- 
tent currents which can be (in the diffusive regime) about 
three-to-four times larger than the white-noise-based result 
Ifyp°^ ~ {evp / L,){1/ L). This finding resembles the exper- 
imental findings of reference [10]. 

Finally, we extend our study to the 3D conductors with 
the columnar grains 01 9142811 . which are fundamentally dif- 
ferent from the tiny randomly-oriented grains, implicitly as- 
sumed in any white-noise-based description of disorder. We 




FIG. 1: Wire made of the 2D conductor of width W and length L. 
The figure on the left depicts the wire with impurities positioned at 
random with random signs of the impurity potentials. The figure on 
the right depicts the impurity-free wire with grain boundaries. The 
grain boundaries are represented by the straight lines which intersect 
the wire with a random orientation in randomly chosen positions. 



show that the typical persistent current in the diffusive metal- 
lic ring with the columnar grains is given by the formula 
Ityp ~ l.2Qs/NH{evF I L){1 / L), where Nh Hkp/Tr is 
the number of the 2D subbands within the thickness H. For 
the Au ring with H = 70nm the formula gives the result 
Ityp ~ 20{evF/ L){1/ L), which is not far from the experi- 
mental results of reference |8]. 

In section II we discuss our calculation of the Landauer 
conductance: we review the scattering-matrix method for the 
wire with impurity disorder and we include the grain bound- 
aries. In section III we describe our calculation of the per- 
sistent current. Our results are discussed in section IV, which 
also contains extension to the case d = 3. The appendices A, 
B, and C describe a few technical aspects. 



II. SCATTERING MATRIX AND CONDUCTANCE 

We consider the electron gas confined in the two- 
dimensional (2D) conductor depicted in the figure [T] At zero 
temperature, the wave function (p{x,y) of the electron at the 
Fermi level (Ep) is described by the Schrodinger equation 



(2) 



with Hamiltonian 



+ + (3) 

where m is the electron effective mass, Ud {x, y) is the poten- 
tial due to disorder, and V{y) is the confining potential due to 
the edges. The confining potential can be expressed as 

rO, 0<y<W 
[ oo, elsewhere 

For the impurity disorder we use the simplest model potential 



(5) 



where we sum over the random impurity positions \xi , yi] with 
a random sign of the impurity strength 7, as it is shown in 
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the figure [T] Disorder due to the grain boundaries can also 
be modeled by means of (|5]), if the individual (5-barriers in 
equation (|5]l are positioned on the straight lines (grain bound- 
aries in the figure [T]) equidistantly and with a positive sign 
of the constant 7. Details will be given later on. Now it is 
important that both the impurity disorder and grain-boundary 
disorder are represented by a sum of the two-dimensional S- 
functions. This allows us to treat both of them by a very simi- 
lar scattering-matrix technique. We first review the scattering- 
matrix technique for the impurity disorder 

Assume that the disordered wire in figure[T|is connected to 
two ballistic semiinfinite contacts of constant width W. In the 
contacts the electrons obey the Schrodinger equation 



(6) 

where V{y) is the confining potential given by equation (|4]i. 
Solving equation ^ one finds the independent solutions 
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Xn{y), ri = l,2,...oo, 



with the wave vectors A;„ given by equation 
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where e„ is the energy of motion in the ?/-direction and 
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0<y <W 
elsewhere 
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is the wave function in the direction y. The vectors fc„ in dTjl 
are assumed to be positive. The energy en+fi^k'^/2m is called 
the n-th energy channel. The channels with e„ < Ep are 
conducting while the channels with e„ > Ep are evanescent. 

We define A^{x) = a^e^**^"^ and B^{x) = h^e^'^^^, 
where and 6^ are the amplitudes of the waves moving in 
the positive and negative directions of the x axis, respectively. 
The wave function ^p{x, y) in the contacts can be expanded in 
the basis of the eigenstates (|7]i. At the boundary a; = 
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^(0,y) = ^[A+(0)+A; 
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while at the boundary x — L 
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^{L,y)^J2[^n{L) + B-{L}]xn{y). dD 

n=l 

where N is the considered number of channels (ideally — 
00). We define the vectors A^(0) and B^{L) with compo- 
nents A^^^ Ar(0) ^iid B^^^ Ar(i)' respectively, and we 



simphfy the notations A (0) and B (i) as A andB . The 
amplitudes A^ and are related through the matrix equation 
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where 
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is the scattering matrix. Its dimensions are 2N x 2N and 
its elements t, r, t', and r' are the matrices with dimensions 
N X N. Physically, t and t' are the transmission amplitudes of 
the waves A"*" andB~, respectively, while r and r' are the cor- 
responding reflection amplitudes. The matrix t is composed of 
the elements tmn{kn), where tmn{kn) is the probability am- 
plitude for transmission from the channel n to channel m. 

If we know im„(fc„), the conductance g can be obtained 
from the Landauer formula ll29ll . In units 2e^/ h it reads 
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where r„ is the transmission probability through disorder for 
the electron impinging disorder in the n-th conducting chan- 
nel and we sum over all (Nc) conducting channels. To obtain 
t„„i, we need to determine the scattering matrix S. 

Consider two wires 1 and 2, described by the scattering ma- 
trices Si and 82- The matrices are defined as 
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Let 
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(16) 



is the scattering matrix of the wire obtained by connecting 
the wires 1 and 2 in series. The matrix is related to the 
matrices Si and 5*2 through the matrix equations jlSll 



ti2 = t2[I - r[r2]-Hi, 

ri2 = ri+t[r2[I -r[r2]-Hi, 

t[2 = t'i[I + r2[I-r[r2]-'r[]t'^, 



(17) 
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t2[I - r[r2]-'r[t'2, 



where I is the unit matrix. The equations ( fTTI l are usually 
written in the symbolic form 



S-_ 
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Si 



(18) 



Consider the wire with impurity potential Q. Between any 
two neighboring impurities there is a region with zero impu- 
rity potential, say the region Xi^i < x < Xi, where the elec- 
tron moves along the x axis like a free particle. The wire with 
n impurities contains n+1 regions with free electron motion, 
separated by n point-like regions where the scattering takes 
place. As illustrated in figure |2] the scattering matrix S of 
such wire can be obtained by applying the combination law 



S ^ pi (E) si (E) P2 <E) S2 ' 



. S„ E)Pn+l, 



(19) 



where pi is the scattering matrix of free motion in the region 
Xi-i < X < Xi and Si is the scattering matrix of the i-th 
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FIG. 2: Wire with the randomly positioned point-like impurities, 
where an individual impurity is represented by the two-dimensional- 
(5-function-like potential with random sign. The n impurities de- 
scribed by the scattering matrices Si divide the wire into the n + 1 
free regions described by the matrices pi. Also shown are the wave 
amplitudes A* andB*. 



impurity. The symbols (E) mean that the composition law (fTsT l 
is applied in (fT9l l step by step: one first combines the matrices 
pi and si, the resulting matrix is combined with p2, etc. 
The scattering matrix pi can be expressed as 



Pi 



$ 
$ 



(20) 



where is the N x N matrix with zero matrix elements and 
$ is the N X N matrix with matrix elements 



Finally, the scattering matrix 
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(22) 



is composed of the matrices 

t = t' = [K 
r = r' = -[K + iVy^iV, 

where K and F are the N x N matrices with matrix elements 

m7 



(23) 
(24) 
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■Xm(y«)xn(y»). 



(25) 



The scattering matrix method for the grain boundaries is the 
same like for the impurities, because a single grain boundary 
is formally modeled by a set of the point-like impurities (see 
figure |3]l. We start with the grain boundaries oriented per- 
pendicularly to the wire. Disorder due to the perpendicular 
boundaries is modeled by the potential 



(26) 



where 70 is the strength of the perpendicular boundary and 
Xi is its random position along the wire. Obviously, the S- 
matrix of the wire with perpendicular boundaries is given by 



the combination law ( |T9b , where Si are the scattering matri- 
ces of the individual boundaries and the matrices pi describe 
the free electron motion between two neighboring boundaries. 
The potential of the perpendicular boundary at a; = reads 



(27) 



Formally, it is a one-dimensional version of the impurity po- 
tential ^5{x)5{y — yi), for which the matrix is known: it is 
given by equations ( |23] ). ( |24| ). and ( |25] |. Therefore, the scat- 
tering matrix Si for the potential ( l27T i is given by the same 
equations, except that the elements of the matrix F now read 



^ nin — ^2 ^ 



(28) 



The elements of the matrices t, t\ r, and r' can be written as 
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t' 
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where = rwyQ/h^. Since the matrices t, t', r, and r' are 
diagonal, the electron impinging the perpendicular boundary 
in the channel n is reflected back to the same channel. The 
reflection probability for the channel n = 1 is 



Rc 



(31) 



where we use the approximation fci ^ kp, with kp being the 
2D Fermi wave vector. In other words, the 2D electron im- 
pinging the grain boundary perpendicularly is reflected from 
kpio ~kp with the reflection probability Rq coinciding with 
|riip. The equation ( |3T1 ) allows us to describe the grain 
boundary by the parameter Rq which is measurable. 

In real metallic wires the perpendicular grain boundaries 
usually do not exist ll30ll . Indeed, the matrix elements ( |29] | and 
( l30l ) are diagonal. This means that there is no inter-channel 
scattering, i.e., transport through such wire takes place in the 
mutually independent channels. However, a sin gle d isordered 
ID channel is always in the localization regime ll3lll32ll while 
the metallic QID wires usually exhibit diffusive regime 13311 . 
So we consider the grain boundaries with random orientation. 

As shown in the figure [3] a single grain boundary is mod- 
eled by a set of the equidistantly-positioned repulsive point- 
like impurities with the nearest-neighbor distance <^ Xp. 
In this model, the potential of the grain boundaries is given as 



where 



UD{x,y) = Y,B,{x,y), (32) 



Bi{x,y)= jS{x - Xij)5{y - yij), (33) 
yj<w 
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FIG. 3; Wire with grain boundaries represented by straight lines. The 
ith grain boundary starts at point qi and ends at point pi . The posi- 
tions Qi and Pi are chosen as random. The avoid appearance of the 
mutually intersecting boundaries, the random numbers are ordered 
increasingly, i.e., qi < q2 ■ ■ ■ < qs and pi < P2 ■ ■ ■ < Ps- The angle 
between the i-th boundary and wire edge is ai, another important 
parameter is the mean lateral size of the grain, do- Inset shows in 
detail a single boundary. The boundary is represented by a set of the 
equidistantly-positioned repulsive point-like impurities (plus signs), 
where a single impurity is modeled as a two-dimensional (5-function- 
like energy barrier (see the text). If we choose the nearest-neighbor 
distance Ac <^ Af, the grain boundary effectively behaves as a 
structure-less one-dimensional energy barrier. 



is the potential of the ith boundary, the same positive 7 is 
used for all impurities, and [xij,yij] is the position of the jth 
impurity at the ith boundary. Following the figure [3] we find 

Xtj =qi+ jAo cos(ai), yij = JAq sm{a.i). (34) 

The grain boundaries described by the potential (I32]l-(l34ll are 
formally a special case of the impurity disorder and therefore 
can be treated by the same scattering-matrix algorithm. 

In our model, the reflectivity of a single randomly-oriented 
grain boundary depends on the parameters Ag and 7. We can 
find the relation between these parameters and parameter Rq, 
defined by equation ( [3TT i. Assume that the grain boundary 
described by potential ( [33] ) intersects the wire perpendicularly 
atx = 0. This simplifies ( [33] ) into the form 



B{x, y) = ^ lHx)d{y - j Ag). 



(35) 



Here jS{x)d{y — jAg) is the same single-impurity poten- 
tial, for which we have akeady expressed the scattering matrix 
(equations |24l|23| andlZSt. Therefore, the matrix F of the po- 
tential ([35] l is simply a sum of the F matrices of all individual 
potentials 6{x)S{y — jAo), i.e., 

(36) 

For small Ag the sum in the equation ([36] ) can be replaced by 
integral X*„i{y)Xn{y)dy and we obtain 

r - ^ ^ X 
Ag 



(37) 



Comparing this expression with (l28l l we obtain the relation 
IG = ll Ag- The perpendicular reflectivity ([3TI 1 becomes 

— , o A 9 — — (38) 




FIG. 4: Circular ring of width W and length L S> W , pierced by 
magnetic flux 

entering the scattering-matrix equation 



Also shown are the wave amplitudes A* and 



where 7 — m^jY?. The randomly-oriented grain boundaries 
can thus be characterized by a single parameter Rq, related to 
the model parameters 7 and Ag through the equation ([38] ). If 
we use Ag <C Af, the resulting wire conductance (for a fixed 
value of Rq) is independent on the choice of 7 and Ag. 



III. CALCULATION OF PERSISTENT CURRENT 

We consider a circular ring of width W and length L ^ W , 
shown in the figure |4| The opening of the ring is pierced by 
magnetic flux <^ due to the magnetic field directed along the 
axis z. The ring is in fact the QID wire from the previous 
text, but circularly shaped and with the wire ends connected. 
Therefore, the electron wave function ip{x,y) and electron en- 
ergy E in the ring can still be described by the Schrodinger 
equation (O with Hamiltonian (O, but we also need to ensure 
the continuity of the wave function and its first derivative at 
the connection. This implies the boundary conditions 



fn{0,y) = exp (-z27r^) (pn{L,y), 



S^i0,y)^e^p[-^2n±)^{^y). 



(39) 



where 0o = h/eis the flux quantum and the exponential fac- 
tor is the Peierls phase factor due to the flux (j). Due to the 
boundary conditions ([39] l the energy E is discrete and depends 
on (f). Now we show how to find the spectrum £'„(</)) 1 14]. 

Since we describe the ring by equations (O and (O, we can 
directly apply the scattering-matrix method developed in the 
preceding text. Indeed, the wave function 1^9(2;, y) can be ex- 
pressed in the ring positions a; = and a; = L by means of 
the expansions ( [T0] | and ([TTT i. where the amplitudes and 
are related through the scattering-matrix equation ([T2] |. If 
we set the expansions ([T0] | and ([TTl i into the boundary condi- 
tions ([39]l. we can rewrite (3% into the matrix form 



k^Al + j^' 



g-i(0) 

Q(0) 



A+ 
B 



(40) 



where Qap{(f>) = exp{i2TT(f)/(f)o)5ap. Combining the matrix 
equations ( fT2b and ( |40] | we obtain the equations 



Q-i(0) 



r t' 
t r' 



which can be rearranged into the form 



t~Q{<t>) 



t'-Q-H<l>) 



0. 



(41) 



(42) 



We label the matrix on the left side of (|42ll as M {E, (/>). To ful- 
fill the equation ( |42] |. the determinant of the matrix M{E, (p) 
has to be zero, i.e.. 



det{M{E,(t>)) =det 



riE) t'{E)^Q-\^) 
t{E)-Q{cP) r'{E) 



= 0. 
(43) 

The submatrices t, r, t' and r' are functions of the electron 
energy E. Therefore, the matrix M is a function of the mag- 
netic flux (f) and energy E. The determinant of M{E, is 
a complex number Therefore, the real as well as imaginary 
parts of det(A/(£', ip)) have to be zero to fulfill the equation 
(l43l l. The equation ( |43] | is thus equivalent to the equation 



|det(Af(i?,0))|2 =0, 



(44) 



which we solve numerically. For a given value of magnetic 
flux, the determinant \det{AI {E , ip))\'^ is calculated numeri- 
cally as a function of the energy E which is varied with a 
small energy step AE from zero up to the Fermi energy. In 
the figure|5] a typical numerical result for |det(A/(i?, is 
shown for a small energy window. The eigen-energies i?„((/>) 
are the zero points of |det(M(i?, In the figure|5]the os- 

cillating function |det(Af (£', 0)) ^ show a series of very sharp 
valleys with a zero minimum value, i.e., the positions of these 
minima are the eigen-energies of interest. We repeat this pro- 
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cedure for magnetic flux (f) + A<j), where 
we obtain the eigen-energies £'„((/> + Aip). 

At zero temperature the persistent current is given as ||2[[3tl 



VE„<Ef 



E 



dE„ 



(45) 



where one sums the single-electron currents for all energy 
levels below the Fermi level. Thus, to obtain the persistent 
current in a single sample, one needs to determine all single- 
electron energies i?„(0 + Ac/)) and En{(f>) below the Fermi 
level, to evaluate all single electrons currents 



dE„ 



(0) 



+ A0) -£;„(0) 



(46) 



and to sum them as shows the equation (1451 ). This procedure 
is computationally cost and allows to study only small rings, 
but later on we introduce a trick allowing to study large rings. 
Finally, we note that describing the ring by means of (|2|l, 
and ( l39b we ignore the magnetic field in the ring. Further, 
the approach ignores the effect of the ring curvature, because 
we assume that the x axis circulates along the ring circumfer- 
ence (see figure|4|i. Both approximations hold for L ^ W. 




FIG. 5: Determinant |det(M(£;, <;/>)) p versus E. The presented data 
points are the values of |det(Af(£', (j)))]^^ calculated for equidistant 
energies with a small energy step AE. The minima of the sharp 
valleys are the zero values of |det(i\/(i5, 4>))\^)- sach zero occurs at 
an eigen-energy En{(j>). All data were obtained for the specific ring 
parameters discussed in figure [Tol but the presented dependence is 
(qualitatively) typical for any disordered ring. 



IV. RESULTS 

A. Conductance of wires with grain boundaries and impurities 

In this subsection we present our scattering-matrix results 
for the wire conductance. The wires with disorder due to the 
randomly-oriented grain boundaries are compared with the 
wires with impurity disorder Also included are the wires with 
disorder due to the perpendicular grain boundaries. 

We use the material parameters m — 9.109 x 10^'^^kg and 
Ep = 5.6eV {Xp — 0.52nm), typical of the Au wires. We 
first study the Au wires of width W = 9nm, with the number 
of the conducting channels being Nc = 34. This number well 
emulates the limit Nc ^ 1, but later we also use larger Nc- 

The parameters of the grain-boundary disorder are chosen 
as follows. The perpendicular reflectivity Rg (equation [38]) 
and the mean lateral size of the grain, da, are kept the same for 
the randomly-oriented as well as perpendicular boundaries in 
order to isolate the effect of random orientation. We recall that 
a single grain boundary is modeled as a line with equidistant 
impurities of strength 7 and nearest-neighbor distance A^. 
Since we keep Aq <^ Xp, the choice of 7 and A^ has no 
effect on the resulting conductance for a given Rq. However, 
once the parameters do, Ac and 7 are chosen, we keep the 
same 7 and the same total number of impurities also in the 
wire with impurity disorder Both types of disorder are thus 
represented by the same numerical model. Therefore, any dif- 
ference between their transport properties reflects exclusively 
the difference between the scattering by repulsive lines and 
scattering by a random array of point-like scatterers. 

In the ensemble of macroscopically identical wires disor- 
der fluctuates from wire to wire and so does the conductance. 
Hence we evaluate (fT4l) for the ensemble of (typically) 10'^ 
wires and we obtain the mean conductance (g), mean resis- 
tance (p) where p — 1/g, and variance (g^) — (g)^. 
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FIG. 6: Mean resistance (p), conductivity o, and conductance fluc- 
tuations ^\&r{g) versus L/l. The results for the wires with impu- 
rity disorder, randomly-oriented grain boundaries, and perpendicular 
grain boundaries are shown by full lines, squares, and dashed lines, 
respectively. The dotted line in the top left panel is the linear fit 
(p) — 1/Nc + pdifL/W, where pdif = 2/(^^0 is the diffusive 
resistivity and the mean free path / is fitted. The conductivity cr is ex- 
tracted from the mean conductance {g) by means of l |48l ) and normal- 
ized by Udi / = 1/ Pdi f ■ The dotted line in the top right panel shows 
the theoretical formula l l50b . The parameters of the grain-boundary 
disorder are Rq — 0.2 and da = lOnm, the impurity strength is 
7 — 0.303, and the impurity density is nj = 2.18nm^^. All three 
types of disorder are shown schematically in the right bottom panel. 



The figure |6] shows our results for the mean resistance, con- 
ductance, and conductance fluctuations in dependence on the 
wire length. If we compare the wires with randomly-oriented 
grain boundaries (data shown by squares) with the wires with 
impurity disorder (data shown by the full lines), we see that 
the results for both types of wires are in good mutual agree- 
ment and also in accord with what one expects for the white- 
noise-Iike disorder. The following features are worth to stress. 

First we look at the mean resistance. Both the impurity 
disorder and randomly-oriented grain boundaries first show 
the linear diffusive dependence lisll 



ki 



Pdif 



(47) 



where pdif is the diffusive resistivity and rtg = fc|,/27r is the 
2D electron density. Notice that in the former case / — 68nm 
while in the latter case I — 13.7nm only. In other words, the 
point-like scatterers constituting the repulsive lines scatter the 
electrons much more effectively like the point-like scatterers 
of the equivalent strength in a random lattice. 

For L I the full line and squares start to deviate from the 
linear rise (l47T l. The deviation is due to the weak localization 
and eventually due to the strong one, manifested by exponen- 
tial rise of (p) with L. On the other hand, for the wire with 
perpendicular grain boundaries (dashed line) we see the expo- 
nential rise of (p) already for L/l ^ 1, which means that the 
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FIG. 7: Typical conductance (Ingr) in dependence on L and L/^. 
The results for the wires with impurity disorder, randomly-oriented 
grain boundaries, and perpendicular grain boundaries are shown by 
full lines, squares, and dashed lines, respectively. The dotted lines 
show the fit (In (7) = — I//^ which gives the localization lengths ^ 
shown in the figure. All parameters are the same as in the figure |6] 



QID wire is in the localization regime. This is because each 
channel behaves like an independent ID disordered channel. 
The figure|6]also shows the wire conductivity 



1 



l/{g)-l/N,W' 



(48) 



normalized by the diffusive conductivity (Jdif = 1/pdif- In 
absence of localization a/adif — 1 independently on L. In 
fact, we see that a /adi / decreases with L linearly both for the 
impurity disorder and randomly-oriented grain boundaries. 
This linear decrease is in accord with the weak-localization- 
mediated behavior predicted for the white-noise-like disorder 
idsllsl]. Indeed, the theory fisllsl] predicts 



(.9) 



(49) 



where a^i/W/L is the classical diffusion term and the term 
1/3 is the weak localization coiTection typical of the QID 
wire. If we write ( |49] l in terms of the conductivity, we obtain 



= 1 - 



1 kf 



L 



CTdif 



3 TTUpW I 



(50) 



In the figure |6] this equation is compared with the numerical 
data for a / adi f ■ Indeed, the agreement is very good both for 
the impurities and randomly-oriented grain boundaries. On 
the contrary, for the perpendicular grain boundaries we see 
that (J jodi f decreases with L exponentially. In such QID wire 
there is no weak localization, only the strong one. 

Finally, we look at the conductance fluctuations ^yvai{g) = 

V (3^) ~ (ff)^- For the QID wire with the white-noise-Iike 
disorder the theory predicts the universal value 1 35, 36h 



Vvar(g) = 0.365. (51) 

The figure |6] shows that the impurity disorder and randomly- 
oriented grain boundaries exhibit ■\/var(5) in accord with pre- 
diction ( fSTI ). For the perpendicular boundaries we see a quite 
different y/vai{g) as the diffusive regime is absent. 

The figure|7]shows the typical conductance (In g) versus the 
wire length. For all three types of disorder, our numerical data 
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FIG. 8: Transmission probability (Tn) versus L/^ for the channel 
indices n = 1, 2, . . . A'^c, where Nc = 34. For n ordered increas- 
ingly, the resulting curves are ordered decreasingly: the top curve 
shows (r„=i), the bottom one shows (r„=jVj.). The data are pre- 
sented separately for the impurity disorder, randomly-oriented grain 
boundaries, and perpendicular grain boundaries. All presented data 
originate from the same calculation as the data in figure|6] 



approach at large L the dependence (In 5) — —L/^ ^6, 371. 
This is a sign of the localization. Fitting of the numerical 
data provides the values of ^ shown in the figure. We find 
the result ^/Z ~ 0.9Nc for the impurity disorder as well as 
for the randomly-oriented grain boundaries. The result — 
0.9Nc reasonably agrees with the result ^/l = Nc predicted 
for the white-noise-like disorder 113 811 and with the numerical 
studies for impurity disorder lfl6ll . For the perpendicular grain 
boundaries we find the value ^/l ~ 3.7, which differs from 
the ID result £,/Iid = 1 LiSJ- The difference is due to the fact 
that our / is the mean over many channels. 

The figure [8] shows the numerical data for (T„). The the- 
ory based on the white-noise disorder predicts, that the con- 
ducting channels are equivalent ifjll I39I1 in the sense that 
(Ti) = {T2) • • • = (TatJ. In the figure[8] this equivalency is 
reasonably confirmed for the wire with impurity disorder and 
for the wire with randomly-oriented grain boundaries. Nev- 
ertheless, in the latter case the equivalency is not so good as 
in the former one. This can be understood if we look at the 
sketch of the grain boundaries in figure|3] It is obvious that the 
boundaries with the angles a — or a — tt are very unlikely 
because the mutually intersecting boundaries are prohibited. 
Consequently, the probability distribution of a in the interval 
(0, tt) is not homogenous: it has a broad maximum around 
tt/2. Our scattering-matrix approach works also for the inter- 
secting boundaries, but such study is beyond the scope of this 
paper: In such case the angle distribution tends to be homoge- 
nous in the whole interval (0, vr), which improves the channel 
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FIG. 9: Mean resistance (p) and conductance fluctuations y^var(g) 
versus L/S,. The results for the wires with randomly-oriented grain 
boundaries (shown by symbols) are compared for various values of 
the reflectivity Rq and grain size da. The results for the wire with 
impurity disorder (the same as in figure|6l are presented in a full line. 



equivalency. Disorder with non-intersecting grain boundaries, 
studied here, is typical for the so-called bamboo-like wires 
140,-43.1 , with do > W. Finally, for the perpendicular bound- 
aries the channel equivalency is absent due to the localization. 

In the figure |9] the wires with the randomly-oriented 
grain boundaries are studied for various values of the grain- 
boundary reflectivity Rq and grain size do- The experimen- 
tally measured values of Rq in the polycrystalline wires range 
from 0.1 up to 0.8 in dependence on the fabrication condi- 
tions lb 



, 44-4^. For all Rq and dc considered in the figure 
|9] the resulting mean resistance and conductance fluctuations 
are similar to the results for impurity disorder (full lines), al- 
beit a noticeable quantitative differences emerge with increas- 
ing Ra- The data for the conductance fluctuations suggest 
that the conductance fluctuations in realistic samples might be 
correlated with the measured values of Rq. However, no mat- 
ter what is the value of Rq, the mean resistance of the QID 
wire with the randomly-oriented grain boundaries rises with 
L linearly up to L c^i ^, where ^ ~ NJ, even for Rg as large 
as 0.8. To see the standard diffusive regime for the strongly 
reflecting boundaries is rather surprising. 



B. Persistent current in rings witli grain boundaries and 
impurities 

In this subsection the persistent currents are studied numer- 
ically in the rings with randomly-oriented grain boundaries 
and rings with impurity disorder. Our numerical results are 
compared with the theoretical result (equation[T]) valid for the 
diffusive rings with white-noise-like disorder. 

In the ensemble of the macroscopically-identical disordered 
rings the persistent current (|45] | strongly fluctuates from sam- 
ple to sample. To asses a typical size of the current in a single 
sample, one can calculate the typical persistent current 



Hyp 



(52) 



where < • • • > means the ensemble averaging. In fact, the 
persistent current / fluctuates also in a single ring when the 
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FIG. 10: The top panels show the single-electron current /„ ver- 
sus the eigen-energy En, calculated by using the relation l l46t and 
algorithm described in the figure|5] The bottom panels show the per- 
sistent current I{Ef) ~ X^vb <e where we sum the single- 
electron currents from the top panels. The dependence I{Ep) is cen- 
tered around zero mean current, as expected for flux $ = ±0.25ft/e. 
The parameters of the ring with randomly-oriented grain boundaries 
are W = 9nm, RG{2eV) = 0.2, da = lOnm, l{2eV) = 13.9nm, 
and L/l = 9.4. For the ring with impurity disorder W — 9nm, 
ni = 2.5nm"^, 7 = 0.4, l{2eV) = 16.7nm, and L/l = 9.0. 



number of the electrons (the Fermi energy) is varied. It has 
been found in [50] that averaging over the electron number, 
performed for a single configuration of disorder, leads to the 
same results as the averaging over different configurations and 
number of particles at the same time. In this work we average 
over the electron number (over the Fermi energy) in a single 
disordered sample. This helps to reduce the computational 
time, but for our purposes still not sufficiently. Fortunately, 
we will see soon that the computational time can be further 
decreased remarkably, when the typical persistent current is 
studied for magnetic flux $ = ± 0.25/i/e. In what follows 
we use $ = — 0.25/i/e. 

The figure [TO] shows the single-electron current /„ versus 
En and persistent current / X^vb <Ef ^" versus Ep, 
calculated for the ring with impurity disorder and ring with 
randomly-oriented boundaries. In both cases the ring param- 
eters (see the figure caption) are chosen to give roughly the 
same mean-free path I and ratio L/l. In spite of their chaotic 
nature, the data for I{Ef ) are centered symmetrically around 
zero mean, which is in accord with the theoretical result 115 ill 



(/) =0, $ = ±0.25/i/e 



(53) 



and which we have also verified by calculating the mean nu- 
merically. It is not trivial that the numerical data for / = 
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FIG. 1 1 : The first row of panels shows the single-electron currents /„ 
from the preceding figure for the energy window 5E — 0.2eV cen- 
tered around the energy E — 1.3eV. The second row of panels shows 
the current I' (Ep) — "^^e <e <Ep^^' where we sum over the 
energy levels in the window up to a level En = Ep. The third row 
of panels shows the current P{Ef) = I' {Ep) — (/')> where (/') is 
the mean obtained by averaging over all Ep in the window 5E. Un- 
like I' {Ep), the dependence P{Ep) is centered around zero mean. 
Note that the data points are connected by a full line which serves as 
a guide for eye but obscures the discrete character of the data, seen 
in the previous figure. 



SvB <Ef ^" plotted in dependence on Ep are centered sym- 
metrically around zero mean current. We stress that the I{Ep) 
dependence (the cloud of the data points in the figure [TOb 
would become strongly asymmetric when just a single elec- 
tron level is omitted (mistakenly or intentionally) from the 
sum J2vE <Ef is just this symmetry around zero mean, 
which allows us to calculate the typical current by means of a 
very efficient trick. Now we explain the trick in detail. 

The first row of panels in the figure [TT] shows the single- 
electron currents /„ from figure [TO] once again, but only for 
the energy window SE — 0.2eV centered around the en- 
ergy E — 1.3eV. The second row of panels shows the current 
I' — J2e <e KEp^ri, where we sum over the energy lev- 
els in the window 5E from the minimum £"„ up to £"„ — Ep. 
Notice that the data for I'{Ef) are not centered around zero 
mean. The third row of panels in the figure [TT] shows the cur- 
rent /''(Sf) = /'(-^F)-(-f')' where (/') is the mean obtained 
by averaging over all Ep in the window 5E. In other words, 
the dependence I'^{Ep) is just the dependence I'{Ef), but 
centered around zero mean artificially. Now we are ready to 
examine the typical persistent current. 

The figure[T2]shows the typical persistent current Ityp, cal- 
culated in dependence on the Fermi energy Ep and normal- 
ized to the theoretical value I*^^'"' = 1.6{evp/ L){1/ L). The 
open circles show the numerical data for Ityp = \J (/^), 
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FIG. 12: Typical persistent current Ityp versus Fermi energy Ef, 
normalized to the theoretical value Ityp°^ = l.Q{evF / L){1/ L). The 

circles show Ityp = ^/ (/^), where / — $I]v£;„<Bp ^" persis- 
tent current due to all single-electron currents below the Fermi level 
(figure fTOt and (...) means averaging in the interval SEp = 0.2eV 

around Ef- The squares show I^yp — \J (/'^) — (-^')^' where 
I' = <E < Bp ^" is the current obtained by summing solely 

the single-electron currents from the energy window 5E around Ef, 
as discussed in figure [TT] The values Ityp originating directly from 
the data in figurefTTIare labeled by arrows. 



single-electron currents below the Fermi level. Such calcula- 
tion is computationally cost because one has to determine all 
single-electron eigen-energies £"„ below the Fermi level. 
However, the figure [12] also shows the numerical data 

(squares) for the typical current I^y^ — ^ (/'^) — (/')^, 

where /' — J^e <e <Ef ^" current obtained by 

summing solely the single-electron currents from the energy 
window SE around Ee, as discussed in figure [TT] This ap- 
proach works much faster because it is no longer necessary 
to determine all i?„ below the Fermi level. Indeed, one only 
needs to determine all £"„ in a small energy window 6E cen- 
tered around Ep- The value of SE should be much larger 
than the typical inter-level distance, but keeping 6E ^ Ee 
still saves a lot of computational time. The figure [12] shows 
that the data for I^yp reproduce the data for Ityp very well. 

Moreover, it can be seen that both calculations agree quite 
well with the theoretical value = l.eievp/ L){l/L). 

The exception are the data in the wire with grain boundaries 



at small Fermi energies. These data deviate from I^yp"^ due 
to onset of the localization regime at small Fermi energies. 
(A closer inspection also shows why such deviation is not ob- 
served for the impurity disorder The reason is that the mean 
free path decays with the Fermi energy much slower than in 
the case of the grain-boundary disorder). 

In what follows we speak about the typical current Ityp but 
we in fact evaluate 1^^^. We consider the Au rings with ma- 
terial parameters m = 9.109 x 10~^^kg and Ee = 5.6eV. In 
figure [13] the typical persistent currents in rings with impurity 
disorder and rings with randomly-oriented grain boundaries 
are studied with impact on the length dependence. The nu- 
merical data for Ityp (shown by symbols) are obtained for 
various ring parameters and compared with the theoretical 
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FIG. 13: Top panels: Typical persistent current Ityp versus the ring 
length L, with Ityp normalized to Jo = cvf / L and with L normal- 
ized to the mean free path I. The symbols show our numerical results 
obtained for various parameters listed in the figure. Other parame- 
ters are m — 9.109 x 10~^^kg and Ef = 5.6eV. The theoretical 
result Jtyp°'' = l.G{evF / L){1/ L) is shown in a full line. Bottom 
panels: The same data as in the top panels, but with Ityp normalized 

J t up • 



result Ityp""^ — 1.6{evE/L){l/L). For the impurity dis- 
order the numerical data agree with the formula Ityp°^ = 
1.6{evE/L){l/L) very well and for the grain boundaries the 
agreement is also very good for large enough L/l. We con- 
clude that for large enough L / 1 not only the impurity disorder 
but also the randomly-oriented grain boundaries behave like 
the white-noise-like disorder 

However, it can also be seen, that for L/l as large as ^ 10 — 
20 the typical current in the ring with the randomly-oriented 
grain boundaries can exceed Ityp°^ by a factor of three to four, 
when the grain-boundary reflectivity Rq is large and/or the 
grain size do is small. It is remarkable that this happens for 
the ring lengths for which the corresponding wire resistivity is 
in the diffusive regime (see the left panel of the figure[9]l. Of 
course, the factor of three to four is too small to explain the 
huge persistent currents (^ cve/L) measured |8] in a single 
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Au ring of length L/l ~ 100. However, it is large enough 
to resemble the experiment where the measured typical 
currents exceeded the formula Ilyp°^' — {evp / L){1/ L) about 
two-to-three times. 



C. Extension to the 3D conductors with columnar grains 

So far we have studied the polycrystalline wires/rings made 
of the 2D conductor of finite width (figure [U. It is intuitively 
clear that the obtained results are representative also for the 
polycrystalline wires/rings made of the 3D conductor, if the 
grain boundaries in the conductor are randomly oriented in the 
3D space. To extend our numerical study to such 3D systems 
is therefore not meaningful. 

It is however meaningful to extend our study to the 3D 
wires/rings with columnar grains flo'-Ss^. In particular, we 
would like to pay attention (see figure O to the bamboo-like 
3D wires with the columnar grains separated by the planar 
boundaries oriented randomly with respect to the wire side- 
walls. In reality, the bamboo-like wires | 40,-43il with the 
columnar grains can be viewed as an opposite hmit to the 
polycrystalline wires composed of the tiny 3D grains (with 
typical size much smaller that the wire cross-section) oriented 
randomly in the wire volume. The bamboo-like 3D wires in 
the figure[T4]are a reasonable idealization of the real bamboo- 
like wires, and we will see that the diffusive persistent currents 
in the rings made of such wires are remarkably larger that the 
white-noise-based prediction 7*^^°'' ~ [evp / L){1 / L). 

We assume (figure fl4l i that the wire of the width W and 
thickness H is connected to the semi-infinite contacts. The 
wave function of the electron at the Fermi level is described 
by the 3D Schrodinger equation 



- — A + Vy{y) + V,{z) + VD{x,y) 



(54) 

where A = d"^ /dx^ + d'^/dy'^ + d"^ jdz^, Vd {x,y) is the 
potential of disorder due to the columnar grain boundaries, 
and Vy{y) and Vz{z) are the confinement potentials: 



Vyiy) 



0, <y <W ^/^^^0, Q< z < H 
oo, elsewhere ' ^ \ oo, elsewhere 

(55) 

First we solve ( |54| ) in the contacts, where we keep Vd — 
as is customary in the Landauer conductance theory lIlSll . For 
Vd = the energies in the directions y and z, and 
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and the corresponding wave functions are 
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(56) 
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FIG. 14: Model of the bamboo-like 3D wire with columnar grains. 
The columnar grain is a grain shaped as a column parallel with the 
growth direction - in our case with the axis z. In the bamboo-like 
wire shown in the figure, the columnar grains are simply the wire seg- 
ments separated by the planar boundaries (shaded areas) randomly 
oriented with respect to the sidewalls. When viewed from the top, 
such 3D wire looks exactly like the 2D wire with the line-shaped 
grain boundaries (figurefTJ, analyzed up to now. 



The wave function in the contacts can thus be expressed as 



oo oo 



(58) 



where the wave vectors kmn obey the equation 



?)2p 
tjp = 

2m 



^rnn: ^mn — ^rn ' (59) 



with being the bottom energy of the channel [m, n]. 
Clearly, kmn is the Fermi wave vector in the channel [m, n]. 
The vectors kmn are real for Ep > emn and imaginary for 
Ep < emn- The number of the conducting channels (channels 
with emn < Ep) is Nc = ttNwNh/'^, where Nw = kpW/ir 
and Nh — kpH/ir are the numbers of the conducting chan- 
nels in the y and z directions, respectively. We can order the 
terms [mn] in the sum dSSl l so that the energies emn are or- 
dered increasingly starting by en. Then the first Nc terms in 
the sum (fSSl ) are due to the conducting channels. 

Now we analyze (l54b in the disordered region, where 
VD{x,y) is not zero. Since VD{x,y) is z-independent, 
it is useful to set into ( |54| | the expansion tp{r) — 
En'=i ^n'ix, y)Xn'iz)- Multiplying (El by x" i^)* and in- 
tegrating over z we get the equations 
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(57) 



Vn{x,y) = 

= Enipn{x,y), (60) 

where n ~ 1,2,... and En — Ep — e^ are the Fermi ener- 
gies in the 2D subbands arising in the vertical direction. The 
equation ( l54b thus splits into a set of equations ( l60l l which are 
formally the same as the Schrodinger equation (|2]i for the 2D 
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conductor. Hence, the disordered 3D conductor in figure [14] 
can be viewed as a parallel connection of independent 2D con- 
ductors with the same disordered potential Vu {x, y), but with 
various Fermi energies £"„. Of course, these 2D conductors 
are in fact the 2D subbands in the vertical direction. 

Therefore, the Landauer conductance of the wire with 
columnar grains can be expressed as 



N„ 



where 



9n{En) 



km {En) 



(61) 



(62) 



is the Landauer conductance ( fT4b rewritten for the n-th 2D 
conductor (?i-th vertical 2D subband) with Fermi energy En. 
We recall that kj{En) — kjn, where kjn is the Fermi wave 
vector in the ID channel [j, n], defined by equations ( fSSl l and 
i59[ . The transmission amplitudes tjm, describing the elec- 
tron transmission through the columnar grain boundaries, can 
be evaluated by means of the same scattering matrix as we 
have introduced in section 11 for the line-shaped grain bound- 
aries (figure O, except that now the Fermi energy is £"„. 

In the figure [15] transport in the 3D wire with columnar 
grains is compared with transport in the corresponding 2D 
wire, obtained from the 3D wire by setting H ^ and 
Nh = 1, and by keeping the same Fermi energy. This means 
that the grain boundaries in the 2D wire are the randomly- 
oriented line-shaped boundaries studied in the preceding text 
(right sketch in the figure [T]). The comparative study shows a 
few results which are worth to stress. 

First a comment on the localization length ^ in the figure 
[TS] Note that the values of ^ in the 3D wire and 2D wire are 
the same. The 3D wire with the columnar grains is a parallel 
connection of independent 2D wires (2D subbands) with 
the same disorder and different Fermi energies. The 2D wire 
with the largest Fermi energy provides the largest localization 
length and this is just the localization length of the whole 3D 
wire because the conductance contributions from other Nh — 1 
wires become negligible for large L. The 2D wire with the 
largest Fermi energy is just the 2D wire obtained from the 3D 
wire by setting — > and by keeping the same Fermi energy. 
As a result, ^ is the same in the 3D and 2D wires. 

Notice now the mean resistance for L < ^. It is roughly 
Nh times smaller for the 3D wire than for the 2D wire and 
the resulting mean-free paths I3D and I2D give the numerical 
ratio l^D /I2D — 0.9, which is in good accord with the formula 



3D 



9tt 
32 



2D 



0.88/^ 



2D, 



(63) 



derived in the appendix B. One also sees that for L < the 
conductivity cr exhibits in the 3D wire essentially the same 
weak localization behavior as in the 2D wire. In summary, 
the resistance and conductance of the 3D wire with columnar 



3D wire with columnar 
grains (Fig. 14) 



2D wire with line-shaped 
grain boundaries (Fig.l) 
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FIG. 15: Mean resistance (p), conductivity cr, and conductance fluc- 
tuations \/ var((?) versus L/^ for the Au wire with grain boundaries. 
The left column of panels shows the results for the 3D wire with the 
columnar grains (figure O, calculated for three different sets of the 
parameters W , H, dc, and Rq- The right column of panels shows 
the results for the 2D wire with the line-shaped grain boundaries, 
which is "fabricated" from the 3D wire with the columnar grains by 
setting H ^ and Nh = 1, and by keeping the same Fermi energy 
Ef — 5.6eV. In other words, the 2D wire and the first vertical 2D 
subband of the 3D wire are the microscopically-identical 2D con- 
ductors. The dashed lines show the linear fit of the (p) versus L 
dependence, from which we obtain the diffusive resistivity pdi / . The 
mean free paths I^d and I2D are extracted from pdi / and from the 3D 
and 2D Drude-resistivity expressions. The conductivity a is obtained 
from the mean conductance (g) by means of | |48I > and normalized by 



grains exhibit a standard diffusive behavior, similarly as for 
disorder which is white-noise-like in the 3D space. The fact 
that the wire is effectively composed of the Nh independent 
2D wires is reflected by the conductance fluctuations: the fig- 

/Nh times larger than 



ure[T5]shows that y/vm-{g) is roughly 
the standard value, which one expects. 

Consider now the ring made of the 3D wire with the colum- 
nar grains. The ring is composed of the Nh independent 2D 
rings. If the n-th 2D ring carries the persistent current In, 
the total persistent current Icg in th^ ring with the columnar 
grains reads 



Nh 

Icg = / ^ In{En). 



(64) 
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To calculate I„ numerically, we evaluate for each individual 
2D ring the spectrum of all single-electron currents below the 
Fermi level (in the same way as in the figure [TO]l and we sum 
these currents to obtain I„ . After that we evaluate the sum 
(|64] | and we eventually perform averaging to obtain the typical 
cuii-ent = ^(/^^). 

We also estimate analytically. The simplest estimate, 
I^yp ~ \/ N h{ev F / L){l2D I L), assumes that each of the Nh 
rings supports the same typical current (the value of which 
is {evp / L){12d/ L) because the columnar grains create the 
white-noise-like 2D disorder in the plane perpendicular to the 
columns). A more precise estimate (appendix C) gives 

Ig'; = ^ l.2Q^{evFlL)(hDlL), 

(65) 

where 

r^^^ = l.\{evp/L){hD/L) (66) 

is the expression ([T]) for d = 3, the 3D mean free path Ij,d — 
0.88^2D (equation|63ll, and the 2D mean free path is assumed 
in the form I2D — CEp (the constant C is determined by 
fitting the numerically calculated I2D as shows the inset to the 
figure[T6]l. 

In the figure [16] the typical current in the ring with the 
columnar grains is calculated as a function of the Fermi en- 
ergy. It can be seen that the formula i65[ (full line) agrees 
quite well with the numerically calculated (open circles), 
while the formula ( |66] ) (dotted line) underestimates the nu- 
merical data about \^Nh times. Obviously, the formula (|66] | 
holds only for disorder which is white-noise-like in the 3D 
space, not the case for the columnar grains. 

The ring considered in the figure [16] is rather small. For 
the Au ring with H = 70nm the formula ( 1651 ) gives the re- 
sult/typ ~ 20{ev f/L){13d/L). This result resembles the ex- 
periment fs*!, where the diffusive persistent currents in three 
individual Au rings of length L/l ^ 100 exceeded the value 
[evp/L) (/3D / L) one-to-two order of magnitude. 

V. SUMMARY AND CONCLUDING REMARKS 
A. Summary of results 

We have studied mesoscopic transport in the QID wires 
and rings made of a 2D conductor of width W and length 
L ^ W. We have compared transport in an impurity-free 
conductor with grain boundaries with transport in a grain-free 
conductor with impurity disorder. 

The transmission through the disordered conductors was 
calculated by the scattering-matrix method, and the Landauer 
conductance has been obtained. We have also calculated the 
persistent current in the rings threaded by magnetic flux: we 
have incorporated into the scattering-matrix method the flux- 
dependent cyclic boundary conditions and we have introduced 




Ep [eV] 

FIG. 16: Typical persistent current in the 3D ring with the columnar 
grains, lu/p, versus the Fermi energy. The ring dimensions are W = 
9nm, H — 7nm, and L — 130nm (L ^ / for all considered Ep), 
the parameters of the grain boundaries are Rg{Ef = 2eV) — 0.2 
and da = lOnm. The circles show the numerical data, the full line 
represents the estimate i65\ , and the dotted line is the white-noise- 
based result l l66t . Inset shows the mean free path in the 2D wire, ho, 
as a function of the Fermi energy: the squares show the numerical 
data and the dotted line is the linear fit ho ~ CEp- 



a trick allowing to study the persistent currents in rings of al- 
most realistic size (the typical persistent current for magnetic 
flux ±0.25/i/e was extracted from the single-electron ener- 
gies in a narrow window around the Fermi energy). We have 
mainly studied the conductance and persistent current in the 
diffusive transport regime. Our results are the following. 

If the grain boundaries are weakly reflecting, the systems 
with the randomly-oriented grain boundaries show the same 
(mean) conductance and the same (typical) persistent current 
as the systems with impurities. The obtained results also agree 
with the single-particle theories of diffusive transport ||a [tI 
[HI], treating disorder as a white-noise-like potential. 

However, if the grain boundaries are strongly reflecting, the 
rings with the randomly-oriented grain boundaries can exhibit 
in the diffusive regime the typical persistent currents about 
three-to-four times larger than the white-noise-based formula 
Ityp°^ ~ {evp/L){l/L). This finding resembles the exper- 
iment [10], where the typical persistent currents measured 
in the diffusive Au rings were two-to-three times larger than 
I^^-'- c {evp/L){l/L). 

We have also extended our study to the 3D conductors with 
the columnar grains. We have shown that the typical persistent 
current in the diffusive metallic ring with the columnar grains 
is given by the formula Ityp ~ l.2Q^/NH{evp / L){1/ L), 
where Nh — Hkp/ir is the number of the 2D subbands 
within the thickness H. For the Au ring with H = 70nm 
the formula gives the result Ityp ~ 20{evp/ L){1/ L), which 
is not far from the experiment [8], where the diffusive persis- 
tent currents measured in three individual Au rings of length 
L/l ~ 100 were one-to-two orders of magnitude larger than 
IZ"^ ^ {evp/m/L). 



14 



B. Comment on relevance for experiment 

Of course, we cannot conclude that our study is a defi- 
nite explanation of the experiments iQl and lioll . because the 
polycrystalline structure of the experimental samples in these 
works is not known. Moreover, even if we would assume 
that the polycrystalline grains in the experiment Jst] are both 
columnar and strongly reflecting, our study predicts the per- 
sistent current only thirty-to-fifty times larger than the formula 
Ily^ — {evp/ L){1/ L) while the largest experimental value 
in [8] exceeds the value [evp / L){1/ L) almost two hundred 
times. 

Nevertheless, according to our study one should not be sur- 
prised when two experiments 12] confirm the formula 
Ifyp"^ — [ev p / L){1 / L) convincingly and two other exper- 
iments id, \\d^ do not. Our study shows clearly within the 
single-particle picture, that the experimental results can de- 
pend quite strongly on the nature of the polycrystalline grains, 
being very different for different fabrication conditions ca- 
ll- 

The columnar grains are fundamentally different from the 
tiny randomly-oriented grains, implicitly assumed in any 3D 
white-noise-based description of disorder Unlike the tiny 
random grains, the columnar grains produce the white-noise- 
like disorder only in the plane perpendicular to the columns, 
which gives rise to the factor -/Wr in the formula Ityp ~ 
1. 26 ^/Nh {ev p / L) {I / L) but which has essentially no effect 
on the diffusive resistance of the wire. It might be instruc- 
tive to fabricate intentionally the diffusive normal-metal rings 
with various types of grains, to measure the persistent cur- 
rent, and to correlate the data with the grain properties: the re- 
sults Ityp ~ {evp/L){l/L) and Ityp ~ \fN^{evF / L){1/ L) 
should appear for the tiny random grains and columnar grains, 
respectively. 

C. Comment on large persistent currents in rings with 
perpendicular grain boundaries 

Finally, we make a comment on the rings and wires with the 
grain boundaries perpendicular to the current. By considering 
the perpendicular boundaries, the work [52] predicted the per- 
sistent currents of size ^ evp jL in rings of length L ^ li^g^ 
where is the 'mean free path' defined as 

or 

h.,-^^9\ (67) 

with (g) being the corresponding wire conductance. Since the 
result ~ evp/L strongly resembles the experimental results 
H, we revisit it briefly. We apply our 2D model (figure[3]l. 

In the figure [T7b we compare the typical persistent currents 
in the rings with the randomly oriented boundaries (empty 
squares) and rings with the perpendicular boundaries (full 
squares), while in the figure Wlh we show the corresponding 
wire conductances (g) . Note that these data are plotted in de- 



^ 3 ^1 ' \ I \ ' \ ' I 

a) : , b)- 




FIG. 17: Panel a: Typical persistent current Ityp versus the ring 
length L, normalized to Jo = evp/L. The empty squares are the 
numerical data for the randomly-oriented grain boundaries. The full 
squares are the numerical data for the perpendicular grain bound- 
aries. The full line shows the formula 11^^°'' = l.Q{evF / L){1/ L). 
The dotted line shows the formula ( I68l >. where ^ = 3.71 is the lo- 
calization length determined earlier. The mean free path / in these 
formulae originates from the wire resistivity pdif = 2/kFl, de- 
termined in subsection IV.A. Panel b: Inverse mean conductance 
1 / (g) versus the wire length L in the wires corresponding to the 
rings in the panel a. The dashed line is the diffusive dependence 
l/{g) = 1/Nc + {2/kFl){L/W). Panel c: The same numerical data 
as in the panel a, but versus the ratio L/lt^g), where /(g) is extracted 
from the data in panel b by means of the formula ( I67t . 



pendence on the length L. We set the conductances from fig- 
urefTTt) into the formula WT\ and we evaluate Z^^^ and Ljli^gy 
In the figurefT?]: we plot the typical currents from figurefTTb in 
dependence on Ljli^gy Indeed, for the perpendicular bound- 
aries we observe the effect predicted in [52J, namely the per- 
sistent currents ~ evpJL for the ring lengths Iijli^g) ^ 1. 
However, in reference [52] this result was not compared with 
the result for the randomly oriented boundaries. The figure 
[T7b shows that the large current diminishes when the orienta- 
tion of the boundaries becomes random. 

Since the values ^ evpjh resemble the large persistent cur- 
rents in the Au rings of the experiment [8], one might spec- 
ulate about presence of the perpendicular grain boundaries in 
the measured rings. However, we see in the figure fTTb that 
the wire with the perpendicular boundaries is in the localiza- 
tion regime for all wire lengths L, for which we observe the 
currents ^ evpjhm the ring. In contrast to this, the Au wires 
used to determine the mean free path experimentally [SJ were 
safely in the diffusive regime, not in the localization regime. 

Moreover, the formula Wh should not be used when the 
wire is in the localization regime. Indeed, the mean free path 
should be length independent, while Z^g^ depends on L quite 
strongly due to the exponentially raising 1 / (5) . 

The figure [T7b also shows that the numerical data for the 
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perpendicular boundaries roughly agree with the formula 

hyp = ^/Kievp/L) exp(-L/20 (68) 

while the numerical data for the randomly oriented bound- 
aries approach the formula lfyp°'^ — 1.6{evF/ L){1/ L). 



The formula (1681 1 describes the QID ring with Nc non- 
communicating channels in the localization regime |4], which 
is not the case for the QID rings of the experiment i^. 



APPENDIX B: MEAN FREE PATH IN 3D WIRE WITH 
COLUMNAR GRAINS 



We set into the equation ( [6T| i the formulae g 
and gn = ^o-fi^, where a^o is the 3D conductivity and af^ 
is the conductivity of the n-th 2D wire. We obtain the equation 



D 



-I Nh 



(74) 
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APPENDIX A: TYPICAL PERSISTENT CURRENT FOR 
MAGNETIC FLUX $ = ±0.25$o 

The persistent current / in a single mesoscopic QID ring 
is a periodic function of magnetic flux $. The period of the 
function is ^ h/e. Therefore, it can be expanded into the 
Fourier series as 



/($) =^/pSin(27rp$/$o) 



(69) 



The current / in a disordered ring strongly fluctuates from 
sample to sample due to the microscopic fluctuations of dis- 
order We are therefore interested in the typical persistent cur- 
rent Ityp — \J (/^), where (...) means averaging over differ- 
ent configurations of disorder Assuming the white-noise-like 
disorder (see the text), the authors of the work iQl derived the 
equation 

oo 

= sin2(27rp<&/$o), (70) 



where 



We set into ( |74l l the Drude expressions 



kp3D 



2D 



(75) 



where kp and l^jj are the 3D Fermi wave vector and 3D mean 
free path, and fc„ = \/2mE„/fi and In are the Fermi wave 
vector and mean free path in the n-th 2D wire. We find that 



3D 



Nh 

E 



(76) 



For the 2D wire with the grain boundaries we expect the de- 
pendence /„ cx E„ - This dependence is in good accord with 
our numerical calculations (see the inset to the figure [T6]l. We 
therefore set into ( |76] | the formula /„ — CEn, where C is 
a constant. Moreover, we also set into (|76] | the expression 
En/Ep = 1 — /Nfj. We obtain the equation 



ho — 



2Nh 



n=l 



1 - 



3/2 



(77) 



where we have used ln=i — CEi ~ CEp = I2D (here I2D 
is the 2D mean free path in the 2D wire with the same Ep as 
in the 3D wire). For Nh ^ 1 the sum in ( |77] i can be replaced 
by integral and we obtain the result 



/ - ^^1 



0.881 



2D, 



(78) 



which is in good accord with our simulation (figurefTSt. 



{I'p) 



96 

7^2p3 



e 

TD 



(71) 



is the mean square of the p-th harmonics, Tp) — L^/D is the 
electron diffusion time around the ring, and D = vpl/dis the 
electron diffusion coefficient. For $ = ±0.25$o the expres- 
sion dTOl l can be rewritten as 



TD 



1 



2 00 

g(2pTiF- 



(72) 



We perform summation in ( |72] | and we obtain the typical per- 
sistent current in the form 



jtheor 
typ 



TD 



(73) 



If we set into ( l73T l the above mentioned expressions for Tp) 
and D, we obtain the equation ([T]i. 



APPENDIX C: TYPICAL PERSISTENT CURRENT IN 3D 
RING WITH COLUMNAR GRAINS 

Using the formula ( |64] |. the typical persistent current I^p 
in the 3D ring with columnar grains can be written as 



tCG 
typ 



\ 



Nh 



(79) 



where we have utilized the fact that the persistent currents I„ 
in the constituting 2D rings are mutually uncorrelated. Due to 
the columnar grains, each 2D ring is subjected to the white- 
noise-like 2D disorder and therefore carries the typical current 
given by the formula ([T]) with d = 2. Thus 



1 /I ^"/fc ''Tt 



(80) 
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We set the last equation into (|79] l and we also apply the equa- 
tion vfj^ = 2En/m and equation Z„ = C£'„ from the preced- 
ing appendix. We obtain the equation 



tCG 
typ 



i.e^cJ- 



Nh 



\ n=l 



(81) 



which can be rewritten into the form 



tCG _ 1 a'^'^F ho 

hyp -^-^ L i 



Nh 

E 

n=l 



(82) 



by using the same procedure as in the preceding appendix. 
For Nh 3> 1 the sum in ( [82] i can be replaced by integral and 
calculated analytically. We obtain the formula 



/35 • L L 



(83) 



which relates the typical current in the 3D ring with the colum- 
nar grains to the typical current in the 2D ring with the same 
Fermi energy and disorder. By means of the formula (fTSl l one 
can rewrite (|83]| into the form ( 
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